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Quantum states and measurements exhibit 
wave-like — continuous, or particle-like — dis- 
crete, character. Hybrid discrete-continuous pho- 
tonic systems are key to investigating fundamen- 
tal quantum phenomena [1—3], generating super- 
positions of macroscopic states [4], and form es- 
sential resources for quantum-enhanced applica- 
tions [5], e.g. entanglement distillation [6, 7] and 
quantum computation [8], as well as highly ef- 
ficient optical telecommunications [9, 10]. Real- 
izing the full potential of these hybrid systems 
requires quantum-optical measurements sensi- 
tive to complementary observables such as field 
quadrature amplitude and photon number [11— 
13]. However, a thorough understanding of the 
practical performance of an optical detector in- 
terpolating between these two regions is absent. 
Here, we report the implementation of full quan- 
tum detector tomography, enabling the charac- 
terization of the simultaneous wave and photon- 
number sensitivities of quantum-optical detec- 
tors. This yields the largest parametrization 
to-date in quantum tomography experiments, 
requiring the development of novel theoretical 
tools. Our results reveal the role of coherence 
in quantum measurements and demonstrate the 
tunability of hybrid quantum-optical detectors. 

Accurate knowledge of a quantum-optical detector is 
essential for its fruitful utilization, be it in foundational 
investigations or technological applications. Photode- 
tectors are normally characterized by several parame- 
ters, including detectivity, spectral sensitivity and noise- 
equivalent power [14]. For quantum detectors, addi- 
tional information is required for a complete specifica- 
tion of the detector. This information is the set of op- 
erators that link the input quantum state of the light 
field to the classical detector output, known as postive- 
operator- valued measure (POVM). It may be estimated 
by means of quantum detector tomogrpahy (QDT) [15- 
18], and is needed if the detector is to be used reli- 
ably. To date, QDT has been successfully applied to 
avalanche photodiodes (APDs) [19], time-multiplexed de- 
tectors [18, 20, 21], transition-edge sensors [22], and 



super-conducting nanowire detectors [23]. The matrix 
representations of the POVMs for these detectors are di- 
agonal in the photon-number basis. Consequently the 
reconstruction problem is linear and positive, and there- 
fore amenable to solution by means familiar to classical 
signal processing [24] . This is not true for a general quan- 
tum detector: the POVM elements can have non-zero off- 
diagonals due to coherent superpositions. Even in con- 
ventional optical communications, coherent modulation 
and detection can increase the data transmission rate by 
an order of magnitude. Moreover, exploration and uti- 
lization of the full Hilbert space of a quantum system re- 
quires a detector capable of implementing a tomograph- 
ically complete set of measurements [25]. Such a capa- 
bility is also vital to fully harness the potential of hybrid 
quantum systems operating at the confluence of discrete 
and continuous variable regimes. To this end, phase- 
sensitive detectors that can measure coherent superposi- 
tions of photon-number states are essential to both quan- 
tum and classical optical applications. We focus on the 
particularly interesting example of the weak-homodyne 
photon detector, which allows us to navigate the discrete- 
continuous confluence region in a tunable manner. 

In this Letter, we introduce a QDT method for the 
reconstruction of the POVM of a coherent optical de- 
tector. This method is applied to two variants of a 
weak-homodyne detector: photon-counting and photon- 
number-resolving (PNRD). The POVM elements of such 
detectors have both phase and number sensitivity, de- 
noted by their off-diagonal and diagonal matrix ele- 
ments respectively. Our experimental procedure, shown 
schematically in Fig. (1), can be universally applied to 
any optical detector. It uses only classical optical states 
as probes, yet, with the resulting POVM we can predict 
the detector response to any quantum state, including 
non-classical ones. Full quantum tomography is only re- 
alized by the development of a new recursive algorithm 
that radically reduces the computational complexity of 
reconstructing the POVM. The new recipe changes the 
complexity from quadratic to linear per recursion in the 
dimension d of the POVM elements. This enables us to 
reconstruct a matrix of unprecedented size representing a 
quantum operation, essential for situations where the size 
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of the Hilbert space is incompletely known. In particu- 
lar, we reconstruct a POVM with 1.8 x 10 6 parameters, 
almost two orders of larger than the largest quantum to- 
mography ever performed. By defining the transformed 
version of the Husimi distribution (Q function), we cast 
the reconstruction problem as a tractable semi-definite 
program, allowing us to determine both diagonal and 
off-diagonal elements. This enables us to characterize 
the detectors' tunability between registering the particle 
and wave behaviors of an input quantum state, making 
our experiment the first full QDT. This also allows us to 
completely characterize phase-sensitive optical detectors. 




FIG. 1: Experimental setup. The output of a laser goes 
through a phase modulator (PM) and an amplitude modu- 
lator (AM) to prepare a set of probe states which are injected 
into the detector. The magnitudes of the probe states \a\ 
are controlled by a motorized half-wave plate followed by a 
Glan-Thompson polarizer. The phases of the probe states 
8 are set by a piezo translator (see Methods: Experimen- 
tal details). From the measurement statistics, the detector 
POVM elements are reconstructed. HWP: half-wave plate. 
GT: Glan-Thompson polarizer. ND: neutral density filter. 
LO: local oscillator. BS: beam splitter. 

QDT is performed by preparing a set of known probe 
states {pm} incident on a quantum detector and observ- 
ing the detector outcomes. The probability of registering 
outcome n is given by the Born rule 



Pn\ m = Tr(/5 m n n ), 



(1) 



where {II„} is the POVM of the detector. The POVM 
elements can be reconstructed from these equations us- 
ing convex optimization. Here we use coherent states as 
probes since they form a tomographically complete set. 
The set of measured outcome statistics gives the Husimi 
Q function of the detector operator Q n (a) = (a|II n |a)/7r, 
where a = \a\e tS is the complex amplitude of the coher- 
ent probe state. Q n {a) contains all the information about 
the detector, and can be used to make any predictions of 
measurement outcomes. In particular, all detection prob- 
abilities are given by the overlap of this Q function with 
the Glauber-Sudarshan P function of the quantum state 
incident upon the detector [26]. As the P- function of a 
non-classical state is highly singular, typically involving 
derivatives of the Dirac-delta function, it makes the over- 
lap extremely sensitive to noise in Q n (a), rendering its 



use impossible in practice. It is thus necessary to recon- 
struct the POVM elements {Iln} themselves. 

For optical detectors, POVM elements can be written 
in the photon- number basis as n„ = JT k Il£ ,fe |j)(fc|. De- 
tector saturation allows truncation of the Hilbert space at 
a finite number of d— 1 photons, leaving us with d 2 — 1 pa- 
rameters to estimate for each POVM element. In quan- 
tum state tomography, one can reconstruct the full den- 
sity matrix because d is typically small due to the lack of 
bright quantum-optical sources. In contrast, quantum- 
optical detectors can have a large dynamic range, with 
d typically in the range of order 10 2 to 10 5 . This makes 
reconstructing 0(d 2 ) parameters extremely challenging. 
We overcome this problem by using the transformed ver- 
sion of the Husimi distribution 



[ 2 * Q n (a)e- ae d8 = 2e~W 2 V n£ 
Jo 
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This reduced function enables the recursive reconstruc- 
tion of the principle diagonal (I = 0) and then each lead- 
ing off-diagonal (I = 1 ■ ■ ■ d) (See Methods: Reconstruc- 
tion procedure). The number of coefficients calculated 
per POVM element per recursion is now no more than 
d, greatly reducing the complexity. For instance, I = 
describes a phase-averaged coherent state as input, for 
which the detection probabilities involve only the princi- 
ple diagonals of the POVM elements. Although one could 
reconstruct the entire operator using Eq. (2), in many sit- 
uations, losses and phase fluctuations restrict the number 
of significant off-diagonals to I -C d. Moreover, situations 
involving input states with a fixed photon number N, like 
NOON states [27] or Holland-Burnett states [28], require 
only TV leading diagonals of the POVM elements to pre- 
dict all measurement outcomes. 

We apply the above strategy to the tomography of a 
weak-homodyne detector. This detector combines the in- 
put state with a local oscillator (LO), typically a weak 
coherent state |«lo) at a beam splitter (BS) with trans- 
missivity T (65.5% in our experiment), followed by an 
APD or a PNRD. The LO provides a phase reference 
to access coherence within the input states through in- 
terference at the BS. This enables the detector to probe 
both the particle and wave nature of a quantum state by 
tuning \otLo\ [H]- Theoretical models of these detectors 
are affected by the BS ratio, LO amplitude and phase, 
as well as the detailed workings of the APD or PNRD. 
Additional factors such as the degree of mode overlap M. 
between the input states and the LO and different losses 
of the input states and the LO, which are difficult to mea- 
sure accurately, must be incorporated into the theoretical 
model. This makes an empirical characterization of these 
detectors, such as QDT, more rigorous and reliable, thus 
revealing their tunable quantum features. 

We first explore experimentally the role of coherence in 
quantum measurements performed by a weak-homodyne 
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FIG. 2: Experimentally reconstructed 1-click POVM elements of a weak-homodyne APD. Left to right: LO strengths |qlo| 2 = 
5.5,0.8,5.5 photons, and a mode overlap M. — 0.99,0.99,0.16 respectively. Only the relative phase between the input states 
and LO determines the measurement result. We assume the LO phase to be 0, making the matrix elements all real. Top: 
(a)-(c): POVM elements in Fock basis. Decreased effective LO strength M\a^o\ 2 leads to suppressed off-diagonal elements of 
the POVM matrix, which represents the coherence in the detector, clearly demonstrating the ability of the weak-homodyne 
detector to bridge the wave-particle gap through a tunable LO. The off-diagonal elements in (c) are more akin to those in (b) 
than (a) due to the effective LO. The fidelities of the reconstructed POVMs n^ oc with the theoretical predictions II„ , given 

by 7 n = ^Tr-y/ y/l^n^y/u^j /Tr(n t n h )Tr(n; oc ) are over 98%. Bottom: (d)-(f): Wigner functions of the above POVMs. 

The measured displacements T> from the origin are 1.76, 0.56 and 0.64, and the values calculated from Eq. (3) are 1.70, 0.65 
and 0.68. The discrepency between experimental results and theoretical predictions can be explained by limited precisions in 
the measurement of the experimental parameters (T, A4 and |aio|)> as well as the noise in the reconstructed Wigner functions 
which limits the accuracy of estimating T>. 



APD (detection efficiency 39%, measured with a sepa- 
rate QDT) by modulating \a.L,o\ an d -M. The resulting 
variations in the relative magnitudes of the diagonals and 
off-diagonals of the POVM elements govern the number 
and phase sensitivities of the detector, as in Fig. (2 a- 
c). When the detection efficiency is fixed, for perfect 
mode-overlap as in Fig. (2 a,b), reducing the LO strength 
leads to suppressed phase sensitivity. Mode-mismatch 
leads to two incoherent processes: (i) The interference of 
the reduced LO with the input state, the former with a 
strength A^|«lo| 2 - ( u ) Mixing of the remainder of the 
LO, (1 — A / J)|a_Lo| 2 I with vacuum at the input. Pro- 
cess (i) reduces the phase sensitivity, and process (ii) has 



no phase sensitivity. The total response of the detec- 
tor is the convolution of the above two processes. From 
this perspective, we can explain why the off-diagonals 
in Figs. (2 b) and (c) are similar. The different LO 
strength between the two cases is offset by their different 
mode overlaps, resulting in similar effective LO strengths 
■M\ctLo\ 2 - The convolution also leads to a modulation in 
the diagonals. This behavior is present in the different 
diagonals of Figs. (2 b) and (c). 

Our results can be further elucidated through the 
Wigner functions of the POVM elements in Fig. (2 d 
- f). The presence of off-diagonals in an operator leads 
to radial asymmetry in its Wigner function. Since the re- 
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FIG. 3: Reconstructed POVM elements for (a) 1-click and (b) 3-click events of the weak-homodyne time-multiplexed PNRD 
in Fock basis. Again the LO phase is 0. We reconstruct up to the 6 leading off-diagonals, beyond which the coefficients are 
negligible. The corresponding Wigner functions are shown in (c) and (d). The simultaneous phase and number sensitivity is 
apparent from the radial asymmetry of the Wigner functions, as is the non-classicality of the operator, from their negativity. 
The f-click Wigner function has an overlap of 98% with a displaced single-photon state (having experienced 84.3% loss). 



sponse of a detector to an input state is determined by the 
overlap of their respective Wigner functions, the phase 
sensitivity of a POVM element can be inferred immedi- 
ately from the radial asymmetry of its Wigner function. 
At one extreme is the Wigner function corresponding to 
a detection outcome for a standard homodyne detector 
(i.e. one in which the LO is strong and the detectors are 
linear photodiodes), which is a delta plane in the phase 
space. The opposite extreme is the Wigner function of 
a Fock state projection, a radial annulus. The exper- 
imentally reconstructed Wigner functions of our weak- 
homodyne detectors demonstrate the tunability between 
the two extremes, that is, field quadrature and photon 
number measurements. The displacement of the dips of 
the Wigner functions from the origin confirms the phase- 
sensitivity of the POVM elements, as in Figs. (2 d-f). 
This displacement is given by [29] 

v=^m M (3) 

matching our measurement results, which go up to 1.8, 
about 3.6 times the quadrature uncertainties of coher- 
ent states. Futher displacement can be achieved with 
higher LO amplitude. Allied with the phase adjustment, 
this tunable region is sufficient to probe most of the non- 
classical state generated to date [4, 13, 30]. The Wigner 



functions in Figs. (2 d, e) have negative values, which 
clearly demonstrate the non-classicality of these detec- 
tors. In contrast, the Wigner function in Fig. 2 (f) is 
positive everywhere. The loss of the quantum feature is 
due to the increased mode-mismatch, which manifests it- 
self as decoherence due to the two competing incoherent 
processes stated above. These results reiterate the impor- 
tance of QDT when quantum detectors become complex; 
it allows us to account for the external degrees of free- 
dom which are difficult to control and may change the 
detector response in an unexpected way. 

Finally, we apply QDT to a weak-homodyne PNRD. 
Such a detector has been predicted to be very powerful 
for non-classical state preparation and measurement [11]. 
It has been experimentally used in the probe of a non- 
classical Wigner function in a point by point manner [13], 
and a demonstration of sensitivity beyond the standard 
quantum limit of coherent optical communication [10]. 
For the PNRD in our setup, a time-multiplexed detector 
with N = 9 time bins [31, 32], convolution effects, lim- 
ited detection efficiency (24%) and interference with the 
LO require us to extend the Hilbert space to d = 450 
(See Methods: Estimation of d). The total number of 
real parameters involved (N — l)d 2 ~ 1.8 x 10 6 , con- 
siderably exceeds the largest quantum tomography ever 
performed, that of an 8-qubit state with 65536 parame- 
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ters [33]. Our recursive reconstruction method provides 
a tractable solution to this problem which is infeasible 
with the standard approach. Figs. (3 a, b) shows the ex- 
perimentally reconstructed POVM elements of the 1-click 
and 3-click events (both displayed up to d — 80 photons) . 
Their distinctive ranges of sensitivity are evident. The 
Wigner function of the 1-click POVM element, shown 
in Fig. (3 c), has an overlap of 98% with a V = 1.62 
displaced single-photon state (having experienced 84.3% 
loss due to the detection efficiency and interference beam 
splitter) . The negativity of the Wigner function is direct 
evidence of its non-classicality, and professes its suitabil- 
ity for engineering non-classical states. 

Phase-sensitive measurements are crucial to fully ex- 
ploit the fundamental features of quantum physics and 
to optimally utilize optical telecommunications channels. 
Weak-field homodyne detectors with photon-number res- 
olution are a unique phase-sensitive measurement in that 
they respond concurrently to both wave- like and particle- 
like characteristics of input quantum states. They hold 
great potential for applications in quantum information 
science operating in a hybrid continuous-discrete setting, 
and fundamental investigations of quantum mechanics. 
We used QDT to elucidate the simultaneous wave and 
particle sensitivity of wcak-homodyne photon-number- 
resolving detection. Our QDT scheme does not rely on 
the technical details of the measurement process, provid- 
ing a universal or device-independent understanding of 
the role of quantum coherence in a measurement process. 
It foreshadows a new means of assessment and verifica- 
tion of more complex optical detectors, for example, that 
can attain superadditive capacity and the Holevo limit 
in coherent communication [9]. 



METHODS 

Experimental details: The local oscillator and 
probe states are generated by an amplified TkSapphire 
laser (Coherent Mira Seed, followed by Coherent RegA 
regenerative amplifier. Operating wavelength Ao = 830 
nm at repetition rate /r = 256.752 kHz). We use a Sem- 
rock interference filter with full-width at half-maximum 
bandwidth of AA = 3 nm to reduce the bandwidth and 
thus the effects arising from dispersion in the optical el- 
ements. 

The laser beam is split at a broadband beam splitter 
(BBS, reflectivity 35%) to generate LO and probe beams. 
The probe and LO then each undergo different control 
elements, to be finally interfered at another BBS. LO 
control is performed by a half- wave plate (HWP) and 
a polarizing beam splitter (PBS), setting the power of 
the LO, as well as defining its polarization. While the 
LO strength is fixed for an experimental run, the probe 
states' amplitude and phase are varied. We use a HWP 
followed by a Glan-Thompson polarizer (GT) to adjust 



the amplitude \a\ of the probe states with a dynamic 
range of 10 5 . A beam sampler with low reflectivity sends 
a fraction of the probe beam to a NIST-traceable Coher- 
ent FieldMaxII-TO power meter to monitor the variable 
attenuation realized by the HWP-GT combination. 

We control the phase of the probe state 9 with a vari- 
able delay line driven by a piezo translator (Physik In- 
strumente P-841.30). Moving the delay line changes the 
phase between LO and probe at the recombination BBS, 
where LO and probe interfere. 

One output of the interference BBS passes through a 
set of pre-calibrated neutral density (ND) filters, and is 
coupled into a single-mode fiber to be detected by an 
APD or time- multiplexed PNRD. The other output of 
the recombination BBS is sent to a fast photodiode to 
monitor 9. The mode-overlap M is also measured with 
this photodiode by balancing the probe and LO, and cal- 
culating the visibility of the interference fringes when 9 
is scanned. 

For each probe state amplitude |a|, we use forty phase 
settings 9, uniformly distributed between and 27r and 
measure the click statistics for 0.5 second. 

Derivation of Eq. (2): Using a = \a\e te , the coher- 
ent state projector can be expressed as 



r,s 

- e - |Q|2 E 
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H r+S i(r-s)8 
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Since Q n (a) = Tr(li n |a)(a|)/7r, with n n = 
T,j,k U n k \j)( k \, Eq. (2) follows from jjf* e^-'-WdB = 

Reconstruction procedure: For an iV-outcome de- 
tector, and p coherent states as probes (|ai), • • • , |a P )), 
we have Np linear equations provided by Eq. (1) p n \ m — 
Tr(|a m )(a m |n„) with m — 1 ■ ■ ■ p and n = 1 • • • N. When 
the maximum photon number is truncated at d — 1, the 
set of linear equations can be rewritten as 



P = FIT, 



(5) 



where P pX N consists of the measurement statistics, F pX( i2 
has as its rows the coherent state probes with F^j — 
e~\ a ^ 2 \ai\ 2 ^ / j\, and tl d 2 xN has as its columns the POVM 
elements that are to be reconstructed (first column has 
the matrix elements of the 0-click POVM, second col- 
umn for the 1-click POVM, etc.). The physical POVM 
set consistent with the data can be estimated using the 
constrained convex optimization 



rmn{\\P-Ftl\\ 2 + g(Il)}, 

JV-l 

subject to tl n > 0, n„ 



(6) 
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where ||M||2 = -\/Tr(MtM) is the Frobenius or the 
Hilbcrt-Schmidt norm, and g(tl) is the regularization 
function whose form will be given later. 

The reconstruction proceeds recursively, starting with 
the diagonals of Ii n . We construct the matrix P l by av- 
eraging the measured statistics as per the left-hand side 
of Eq. (2). A similar averaging of the input states leads 
to an input matrix F l , and satisfies the identity 

P G) = jF «n« (7) 

where IT') corresponds to only the I th diagonal of the 
POVM matrices. For each I, one can then setup an in- 
dependent semi-definite program as in Eq. (6). The con- 
straints however need to be tailored to this new scenario. 
For I — 0, this is trivial, as the diagonal elements must 
all be positive, and the sum over the POVM elements be 
the identity 

N-l 
n=Q 

This same condition requires that the corresponding off- 
diagonal elements must all add to zero. The positivity 
condition is recursively enforced using Sylvester's crite- 
rion, which states that a matrix is positive if and only if 
all of its principal minors are positive. As an example, 
for I = 1 the condition is 

I : 

|n«' +1 | < Vn^n?+ lj+1 Vn, j; £ n« +1 = o Vj. 

(9) 

Beginning with a positive diagonal matrix, which we ob- 
tain by solving the I = reconstruction, the successive 
off-diagonals are added, ensuring at each step that the 
operator thus obtained is positive. 

The reconstruction problem effectively deconvolves a 
coherent state from the statistics to obtain the POVM 
set. This is an ill-conditioned problem, as seen by the 
large ratio between the largest and smallest singular val- 
ues of the matrix F. This makes the POVM extremely 
vulnerable to small fluctuations in the statistics. This 
instability is taken care of by the regularization function 
g, a convex quadratic function, that still allows us to cast 
the regularized problem as a semi-definite program [20]. 
The same regularization function is enforced for each I, 
one that penalizes large differences |n^' J+z — fl^ + +1 | 
as g(IL) = -yJ2j \^n j+l -& n +1 ' j+l+1 \ 2 - Note that the reg- 
ularization makes no assumption about the details of the 
quantum detector or the actual value of 7. Variations of 
7 over two orders of magnitude produced only around 
10% difference in the reconstructed POVM. 

Estimation of d: The response of the time- 
multiplexed PNRD can be modeled as a loss matrix fol- 
lowed by a convolution matrix, which accounts for the 
effect of the beam-splitter network. As in Ref. [20], it 



can be shown that without losses, saturating 8 bins with 
a probability of 99% requires 52 photons or more (as- 
suming both beam-splitters are 50:50). Limited trans- 
mission (0.157, the transmission of the interference BS 
0.655 times the efficiency of the time-multiplexed PNRD 
0.24) adds other statistical reduction to the number of 
photons. To have at least 52 photons survive the losses 
with a probability of 99% requires more than 438 photons 
at the input. Finally, taking into account the destructive 
interference with LO, it requires approximately 450 pho- 
tons to saturate the weak-homodyne time-multiplexed 
PNRD. 
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